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System for Encoding and Manipulating Models of Objects 
Field of the Invention 

The present invention relates generally to computerized systems that 
generate models from the geometry, shape, and other attributes of objects, 
5 and, more particularly, to systems that manipulate the models. The models 
are usable by computer-aided design tools, robotics, molecular modeling, 
kinematics and dynamics of linked bodies, collision detection, satellite 
positioning, articulated characters in animation and for use in the 
educational, scientific, business, and entertainment fields. 

10 

Background 

Points, lines, planes, and spheres of all dimensions can be used to describe 
basic objects or components. These basic objects can be assembles into more 
15 complex objects. Therefore, the manipulation and proximity testing of the 
basic objects has wide-ranging implications for many problems in 
engineering, computer graphics, and the practical sciences. 

For example, in robotics and many haptic devices, a tool is typically 
20 composed of rigid links connected to each other at joints. The possible 
positions of the distal end of each link, assuming complete rotational 
freedom at the joints, is of course a sphere, see U.S. Patents 6,070,109, 
6,064,168, and 5,973,678. 
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In order to check for interference with an obstacle in the vicinity of the tool, 
one can first check to see if bounding spheres intersects. This is a much 
simpler computation than determining exactly the coordinates of all of the 
individual parts. Only when the bounding spheres intersect does it become 
5 necessary to do the more rigorous computation on the exact coordinates to 
determine actual interference. 

In rigid body kinematics and CAD/CAM, a frequent problem is to check for 
collisions between arbitrarily shaped objects, see U.S. Patents 6,054,997 and 
10 5,943,065. This problem can also be simplified by first circumscribing each 

□ object with a bounding sphere. Now, one can first simply check for 

m collisions between spheres, and only when there is a possible collision of 

3 spheres does one have to do the more complex computation to check for 

□ actual collision. 

r 15 

In graphics rendering, it is common to parametrically represent objects as 
^ polygons or Bezier surface patches. Frequently, it is necessary to rotate 

y objects. The rotation of a complex object can simply be done by rotating 

points on the patches (planes) along spherical surfaces having their origin at 
20 the center of rotation. 



In fluid dynamics, and electrostatic problems, one frequently manipulates 
angle preserving conformal surfaces, see U.S. Patent 5,453,934. In higher 
dimensional geometries, a conformal surface can be represented by spherical 
25 rotation. Clearly, operations with a spherical surface would be much easier 
than operating on a complex surface such an airfoil. 
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In a Cartesian model of Euclidean space, the calculations performed on 
spheres are generally of the form: 

Operator((x+hlf +(y+ilf+(z+jl)-rl 2 , (x+h2) 2 +(y+i2) 2 +(z+j2)-r2 2 ) 
where "Operator" represents some spherical operation such as intersection, 
5 union, containment, tangent, translation, rotation, etc., and combinations 
thereof. In most of these application, event these simpler spherical 
computations, due their enormous numbers, still consume most of the 
system's resources when manipulating object models. Similar calculations 
can be performed on planes. 

10 

□ Therefore, it is desired to provide a new method for representing models of 

m objects which solves these problems for practical modeling and simulation 

3 problems of physical objects. Furthermore, it desired to manipulate the 

h models without having to consider the actual Cartesian coordinates of the 

15 underlying objects. 

::S Summary of the Invention 

The present invention encodes points of components of objects measured in 
20 Euclidean three- space G(3,0) to null vectors in a general homogeneous four- 
space G(4,l). This encoding makes it possible to directly determine scalar 
distances between points by vector inner products. 

In addition, planes and spheres planes in G(3) can be represented by non- 
25 null vectors in G(4J) so that distances between planes can also be computed 
by vector inner products which yields a scalar. 
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All relations among points, lines, planes, and spheres in G(3) can be 
represented or determined in G(4,l) without having to resort to their actual 
coordinates, as in the prior art. 

5 Common operations such as 3-D rotational and translational kinematics of 
objects is described by a single spinor equation in G(4,l). The 3-D rotational 
and translational dynamics of a rigid object is described by a single bivector 
equation in G(4,l). 

10 More particularly, the invention provides a method for modeling an object 
composed of one or more components. Data are input to a memory of a 
computer system for each component of the object. The data of each 
component point include Cartesian coordinates expressed in Euclidean space 
of a plurality of points x of each component. 

15 

Each component point x is encoded as a vector x in a general homogeneous 

1 2 1 2 

space byx = (x+2xe + e*)E = xE-2*e + e*, where e and e* are basis 
null vectors of a Minkowski space E. General homogeneous operators are 
associated with each data point to generate a model of the object. 

20 

The general homogeneous operators are applied to each encoded point of the 
associated component for each component to manipulate the model of the 
object. The operators can be in terms of processes executing in the computer 
system. 
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Brief Description of the Drawings 

Figure 1 is a block diagram of a modeling system according to the invention; 

Figure 2 is a graph of a Minkowski plane used by the present invention; 

Figure 3 is a graph of a homogeneous space where objects are modeled 
according to the invention; 

Figure 4 is a table of transformation operations to be applied to objects 
according to the invention; and 

Figure 5 is a schematic of a robot and an obstacle modeled by the invention. 
Detailed Description of the Preferred Embodiment 
System Overview 

As shown in Figure 1, our invention is a system 100 for generating and 
manipulating a computerized model of a real or virtual object 101, for 
example a robot. The object is composed of basic components such as links 
and joints. The object is the "input" 102 into the system 100. The input can 
be described in terms of the geometry (data or constructs) 103 of the basic 
components, and the operations (methods or operators) 104 that can be 
applied to the basic components. 
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The data 103 can be points, lines, planes, and spheres. The data can be 
acquired by measuring the objects, or the data can be parametrically defined. 
The data are stored in a memory of a computer system. The methods 104 can 
be linkage dynamics, torque forces, transformations, rotations, constraints, 
5 and the like. The methods can be implemented as programs executing in the 
computer. 

The system 100 according to our invention operates the model in a 
homogeneous space (G(4,l)). Homogeneous space is described in greater 

10 detail below. We presume that the methods 104 are already in a 

homogeneous form. However, the data 103 expressing the geometry of the 
basic components most likely use Cartesian coordinates of Euclidean space 
(G(3,0)). Therefore, step 115 can convert the data 103 to homogeneous form 
105 where necessary. The encoder 110 combines the homogeneous data and 

15 corresponding homogeneous methods to generate classes that together for an 
object-oriented programming structure (OOPS) 130. 

The modeler 120 operates on the OOPS 130 using run-time parameters 121 
to determine distances, intersections, and tangencies of the basic 
20 components, and to perform operations such as rotations and displacements 
of the basic components. In an object oriented programming structure, 
methods are applied to data. Output of the modeler 120 can be rendered as 
images on a display device 109. 

25 As will be described, the result of many of the operations, for example, 

operation that measure distances, result in a scalar. Therefore, the result does 
not need to be decoded, and can be directly used in the context of the 
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originating Euclidean space. In other words, as an advantage of our 
invention, the output 108 of the modeler 120 can be directly processed by 
any standard rendering engine, or the robotic tool itself, without decoding or 
converting. 

5 

The homogeneous constructs and operators (OOPS) 130 generated by the 
encoder 110 and used by modeler are derived using geometric algebra. 
Geometric algebra is described extensively by Hestenes in "New 
Foundations for Classical Mechanics" Kluwer (second edition), 1999. 

10 Geometric algebra has been used extensively to describe the mechanics of 
particles, rigid bodies, biomechanics, robotics, computer vision, computer 
aided design, orbital mechanics of space programs, and many other practical 
fields. Geometric algebra is a general modeling system for a wide range of 
physical phenomena. Starting from a small number of simple rules that 

15 define the geometric algebra, one can derive vector algebra, quaternions, and 
spinor representations of quantum mechanics. However, applications of 
geometric algebra in our "generalized" homogeneous space are new. 

The standard model for Euclidean space Ef 1 is an n-dimensional real vector 
20 (Cartesian) space R n or, equivalently, a set of real coordinates. One trouble 
with this model is that the origin is a distinguished element, whereas all the 
otherpoints of £" are identical. This deficiency in the vector space model 
can be corrected by removing the origin from the plane and placing it in a 
space one dimension higher. 

25 
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We do this by placing our model in general homogeneous space. In general 
homogeneous space there are no coordinates. Our model also provides for 
an adequate representation for Euclidean points or lines at infinity. 

The idea underlying our general homogeneous space for "points' 1 in R n is to 
remove the troublesome origin by embedding R n in a space of higher 
dimension with a displacing "null" vector. As stated above, we measure or 
model objects in the Euclidean space E l as a set of points with parametric 
properties. A standard way to do this is by identifying each Euclidean point 
of an object with a vector x in R n , as expressed by the isomorphism 
E n =R\ 

We call this the inhomogeneous model of because its origin 0 is a 
distinguished point in R n , even though all points in E 1 are supposed to be 
identical. 

Constructs 

Following are some of the basic data constructs used by our encoder and 
modeler. 

Scalars and Vectors 

Geometric algebra is built over a real, n-dimensional vector space. Thus, we 
use the structures of scalars and vectors as are well known in the art of 
geometric algebra. 



8 



Merl-1281 
Hestenes et al. 



Null Vector 

If a vector a * 0, but the absolute value of vector \a\ = 0, then vector a is a 
null vector. In geometric algebra a null vector can exist when a ^ 0. 

K-blade 

The outer product of k vectors, determined by 

s defined as the antisymmetric part of the geometric product a } a 2 . . . a k . The 
resulting quantity is called a k-blade, where k is called the grade of the 
blade. 

Pseudoscalar 

The largest grade k-blades are called pseudoscalars. We reserve / as a 
special symbol for a largest k-blades because these are all scalar multiples of 
a single unit pseudoscalar. 

Versor 

A quantity 



U — Ufa Ui<_] . . . u 



(i) 



generated by a product of non-null vectors is called a versor 
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Operators 

In order to properly define a class, it is also necessary to define the 
applicable operators (methods) 104. For the above constructs, the operations 
are defined according to geometric algebra. Thus, for example, a linear 
transformation U on the vector space R M is said to be orthogonal when the 
transformation preserves the inner product of vectors, as specified by: 
(Uay(Ub) = a*b. 

Every orthogonal transformation U is determined by a versor U via the 
equation: 

U(x) = eUxIT 1 , (2) 
where 8 is the parity of U. This reduces orthogonal transformations to 
multiplication by versors. 

A multiplication by versors has at least two great advantages. First, it 
reduces composition of orthogonal transformations to multiplication of 
versors. Second, it enables us to represent and visualize any orthogonal 
transformation as a set of vectors via equation (1). Third, the result is a 
scalar that provides the answer to the problem originally expressed in 
Euclidean space. 

In addition, we also define a new class, which leads to the encoding of the 
input 102 into new data structures 104 that are subsequently operated on 
with new methods. 
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Minkowski plane E 

As shown in Figure 2, a Minkowski E plane 20 lis generated, from two 
unique basis null vectors e 202 and e* 203, i.e., a basis vector {e, e* } with a 

5 fixed direction and variable scale. The null basis vectors have the following 
constraints 

2 r\ 2 

e - 0 = e* , 

and 

ee* = 1 + E, e*e = 1 - £, 

10 

where E = eAe* is a unit 2-blade, that is, E = 1. This is called the 
Minkowski plane £. An alternative but equivalent orthonormal unique basis 
vector, with a variable direction and fixed scale, is {e_,e + } 204. Because the 
orthonormal basis vector 204 has a fixed scale but an arbitrary direction, we 
15 prefer the null basis vectors 202-203. 

Conformal Split 

As shown in Figure 3, a conformal split allows us to relate objects and 
20 operations between spaces of different dimensions, particularly to relate 
operations (methods) on data expressed in Euclidean and homogeneous 
spaces. The conformal split is realized by intersecting the cone of null 
vectors 301 generated by [e, e* } with the plane orthogonal to e, resulting in 
the so-called horosphere 300. The horosphere was originally described by F. 
25 A. Wachter (1792-1817), a student of Gauss. He showed that a certain type 
of surface in hyperbolic geometry is geometrically equivalent to Euclidean 
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space. Our invention shows that the horosphere 300 constitutes a 
homogeneous model of Euclidean space. 

Unlike Wachter, we formulate the horosphere in terms of geometric algebra. 
Thus, our horosphere finally becomes a practical tool many real-world 
modeling applications. 

Homogeneous Point 

Using the horosphere 300, we encode any input data point x of an object 
component in the data 103 in standard Euclidean space, into a new point x in 
a homogeneous space, using 

x = (x + 2* 2 e + ei)E = xE- ix 2 e + e* , (3) 
which expresses x in terms of input x. 

Homogeneous Point, Line, Plane and Hyperplane 

For independent input vectors a 0 , a b . . . a m of the data 103, where m=0 
defines a point, m=l defines a line, m=2 defines a plane, m=3 defines a 
hyperplane, etc., a corresponding vector p is encoded respectively for a 
homogeneous point, line, plane, hyperplane by 
p~ IE (e/\ aoA a]A... Aa m ). 
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Homogeneous Sphere 

An input sphere of a radius r, and a center c having coordinates of m 
dimensions, we encode a m-dimensional vector s, called the homogeneous 
5 sphere, by 

1 2 

s = c + i r e. 

Note that in our new homogeneous classes all objects, i.e., points, lines, 
planes, hyperplanes, and spheres of any dimensions, are exclusively 
10 represented as vectors. 

Thus, our invention leads to much simpler novel operations used by our 
modeler 120. Specifically we describe the following new methods, which 
operate on our homogeneous representation of the basic homogeneous 
is constructs. 

Distance between Points 

The Euclidean distance d ah between two Euclidean input points a and b is 
20 determined for homogeneous encoded points a and b by 

2 2 

d a b = (a-b) = -2a*b, 
a scalar directly useable in Euclidean space without having to do any further 
conversion. 
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The length of line, or line segment through two input points a and b is 
determined for encoded homogeneous points a and b by 
(length) 2 = {eAa,\bf = (a - b)\ 

also a scaler. 
Triangular Area 

The homogeneous area of plane given by input points a, b, c, or the 
triangular plane segment with vertices a, b ? c is determined by 

2 J. 2 

(area) = 4 (eAaAbAc) ) 

a scalar. 

Distances between Points and Objects 

The distances between homogeneous point a, plane p and sphere s are 
determined by the respective inner products 
a*p, a*b, a*s, p m s, 

a scalar. 

Distances between Spheres 



The scalar distances between two homogeneous spheres sj = c } + \ r^e and 
I 2 

^2 = c 2 + 2 r 2 e is determined by 

W2 = O'Q + 2 (>V 2 + r 2 2 ) = -i [(c 7 - c 2 ) 2 - (rj 2 + r 2 2 )]. 



14 



* 



Merl-1281 
Hestenes et al. 



Rigid Body Motion 



According to Charles' Theorem, any rigid displacement can be expressed as 
a screw displacement. Therefore, the motion of a rigid body in homogeneous 
space is determined by a time dependent displacement versor D = D(t) 
satisfying the differential equation 

b =im 

with a "screw velocity" V given by 
y=_/oo + e\, 

where co is the angular velocity and v is the translational velocity of the 
body. 



Dynamics of a Rigid Body 



The dynamics of a rigid body in homogeneous space is determined by the 
differential equation 

P = W 

with 

P = -IL + e*p and W=-/T + **F, 
where L is the angular momentum and p is the translational momentum of 
the body, while T is the net torque and F is the net force on the body. 



Conformal Transformations 



When equation (3) is substituted into (2) we get 

eU(x + ix 2 e + e*)ElT J E = o(x' + W 2 e + e*) 
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where a is a scale factor. The scale factor can be used for normalization 
purposes. This equation determines a transformation of each input vector x 
into a unique vector 
x' = g(x). 

This induced function g(x) is a conformal transformation of input vectors 
and every conformal transformation can be generated in this way. A 
transformation is said to be conformal if it preserves angles between 
intersecting curves in the space being transformed. 

Every conformal transformation function g(x) can be represented by a versor 
U and composition of conformal transformations is reduced to multiplication 
of versors as expressed by equation (2). 

The ability to compute conformal transformations by multiplication this way 
is new and important innovation for manipulation of objects. 

Figure 4 shows all conformal transformations that can be applied to input 
points and homogeneous points. In our homogeneous form, composition of 
transformations is accomplished by simple linear multiplications, as opposed 
to nesting functions in the standard Euclidean form. Also, note that our 
homogeneous transformations never require division or addition. 

In Figure 4, column 401 shows different types of transformations, column 
402 non-linear transformations of prior art Euclidean space, column 403 
linear homogeneous transformations, and column 404 scalars. In Figure 4, 
the symbols are used as follows, x a point in Euclidean space, n normal to a 
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plane, e a null vector, S a displacement, p a radius, c a center, R a rotation, T 
a translation, o and A,. 

Application Example 

Figure 5 shows a robot 501 operating within the confines of an obstacle 502. 
The obstacle is associated with bounding planes and box 503, and bounding 
spheres 504, and the robot is associated with bounding spheres 505. As the 
robot 501 moves, it is important to avoid hitting the obstacle 502. The 
calculation to determine interference must be made quickly. 

Our invention models the robot and obstacle by a few well chosen spheres 
and planes in homogeneous form. We then test for interference using the 
homogeneous operations as defined above. If the spheres that model the 
robot are not allowed to intersect the spheres and planes of the obstacle, then 
it is ensured that the robot will not collide with the obstacle. Therefore, we 
encode the robot and its operation in homogeneous form, and then apply the 
appropriate operations (methods) 104 to determine the distances between the 
various spheres, planes, and lines. 

In addition, our operations can determine the dynamics of each robot link, 
for given input torques, forces and momentums, and link lengths. Similar 
methods apply, for instance, in CAD/CAM where fast intersection of objects 
can use a collection of bounding spheres. 

Finally, the motion (transformations) of the robot can be determined because 
all rigid body motions are conformal. The modeler 120 can determine the 
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motion by a linear sequence of multiplications, instead of the more 
traditional combinations of adds and multiplies. 

Conclusions and Effect of the Invention 

We have described how an object, measured or defined in Euclidean space, 
can be encoded as a model in a homogeneous space. After encoding, 
complex non-linear operations in Euclidean space become simple linear 
operations in homogeneous space.. This greatly simplifies operations on 
planes and spheres that are frequently used in robotics, kinematics, and 
computer graphics. 

Although the invention has been described by way of examples of preferred 
embodiments, it is to be understood that various other adaptations and 
modifications may be made within the spirit and scope of the invention. 
Therefore, it is the object of the appended claims to cover all such variations 
and modifications as come within the true spirit and scope of the invention. 



18 



